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Abstract The Madelung equations map the non-relativistic time-dependent 
Schrodinger equation into hydrodynamic equations of a virtual fluid. Here we show 
that an increase of the Boltzmann entropy of this Madelung fluid is proportional 
to the expectation value of its velocity divergence. Hence, entropy growth is ac¬ 
companied by expansion resulting from the ability of the Madelung fluid to be 
compressible. The compressibility itself reflects superposition of solutions of the 
Schrodinger equation. Thus, in unitary processes where the Madelung fluid ex¬ 
pands and then shrinks, the Boltzmann entropy may, correspondingly, grow and 
then decrease. The notion of entropy growth due to expansion is common in dif¬ 
fusive processes, however in the latter the process is irreversible. Much unlike 
the Boltzmann entropy, the von Neumann entropy, does not vary with time. To 
elucidate the physical underpinning of the Boltzmann entropy, we examine sev¬ 
eral specific examples. We demonstrate that, for classical diffusive processes, the 
“force” accelerating diffusion has the form of the positive gradient of the quantum 
Bohm potential. In the Madelung fluid, the advective and the diffusive veloci¬ 
ties correspond respectively to the the real and imaginary parts of the complex 
momentum. We find that the diffusion coefficient provides a lower bound of Heisen¬ 
berg uncertainty type product between the gas mean free path and the Brownian 
momentum. 
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1 Introduction 


The Madelung equations [T] transform the non-relativistic time-dependent 
Schrddinger equation into hydrodynamical equations of an Eulerian fluid pisllil 
E] . The aim of this study is to explore the entropy properties of the Madelung fluid 
and compare it with the ones of classical Eulerian fluids. 

The hydrodynamical transformation is obtained when considering the Schrddinger 
equation, 


Q\J/ 

ih^ = H^ = 
dt 






( 1 ) 


for a continuous wave function t) = y^(r, (so that p = of 

a particle with mass m, in the presence of an external potential U{r,t). Using 
the de Broglie guiding equation, Ua = V5 (where the tilde superscripts represent 
hereinafter quantities per unit mass m), the real part of (1) becomes the continuity 
equation 

dt 


= -V • (pUa). 


( 2 ) 


The imaginary part of the Schrddinger equation becomes the Eulerian fluid mo¬ 
mentum equation, 

= -VQ{p) - VU. (3) 


Here, ^ ^ + Ua • V, is the material (advective) time derivative of a fluid 

element along its trajectory and Q = ~ 2 ^ Bohm potential per unit 

mass [S]. 

In Ref. [7], it was suggested that the conservation of the domain integrated 
energy in (1) implies that the Madelung fluid is adiabatic. The fluid conserves the 
sum of the domain integrated kinetic, potential and internal energy, where the 
latter is given by the Fisher information. The domain averaged adiabaticy is in 
agreement with the conservation of the von Neumann entropy. The von Neumann 
entropy, i.e., the trace EntvN = —fcs ^’’'[plnp], where ks is the Boltzmann 
constant and p the density matrix of a closed system with a general (possibly time 
dependent) Hamiltonian, cannot (by virtue of unitary time evolution) change with 
time. One may, nevertheless, devise other illuminating entropy functionals (e.g., 
the “diagonal entropy” of El) that, even for closed systems, transparently adhere to 
standard thermodynamic relations (including the second law of thermodynamics). 
In the current work, our focus is on quantities associated with the Madelung fluid. 
Domain averaged adiabaticy of the Madelung fluid does not imply that entropy is 
materially conserved by a “fluid parcel”. This viable non-conservation differs from 
that of the thermodynamic entropy, p = In (T'^”/p^^®), where Cv is the specihc 
heat of an N particle adiabatic ideal gas, and %^fj = 0 (see, e.g., [10]). We will 
hnd that the classical Madelung fluid dynamics indeed motivates the introduction 
of a “Boltzmann entropy of the Madelung fluid” {Ents) and illustrates that the 
dynamics is intuitively appealing. Specihcally, for the standard density p = E*E 
in, e.g., real space, we set 


Ents = —ks j plnpdU = ks J sdV. 


(4) 
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For Joule free expansion of an ideal gas into a vacunm the temperatnre remains 
constant, the standard increase of the thermodynamic entropy is captnred by the 
change in Ents- The non-trivial temporal changes of the Boltzmann entropy of 
(4) are, generally, very different from those of the von Nenmann entropy. The 
Boltzmann entropy is intimately linked to the Shannon information la¬ 
in Section 2, we explicitly discnss the conservation of the von Nenmann entropy 
for a diagonal density matrix p. We nnderscore that by virtue of the invariance of 
the trace defining Svn under nnitary time evolntion, the von Neumann entropy is 
time independent for any initial density matrix (whether diagonal or not). The new 
form of Svn that we derive will lucidly relate it to the Madelung fluid density p. 
In Section 3, we will examine the Boltzmann entropy production in the Madelung 
fluid via expansion by considering examples of both reversible and irreversible 
processes. 

The continuity equation (2) represents advective rather than diffusive dynam¬ 
ics, where the negative value of the Bohm potential gradient acts as a force in the 
Madelung fluid (3). Nevertheless, in Section 4 we show that for classical diffusion 
the positive value of the Bohm potential gradient acts as an effective force. Besides 
the realization of the quantum Bohm potential in a classical phenomenon this im¬ 
plies that the Bohm potential represents diffusive processes in the Madelung fluid, 
as discussed in Section 5. Summary and conclusions appear in Section 6. 


2 The von Neumann entropy conservation in the Madelung fluid 

In the sections that follow we will largely focus on the Boltzmann entropy as¬ 
sociated with the Madelung fluid. Before doing so, we will now briefly derive 
the von Neumann entropy Svn for this fluid. For a continuous wave function 
I'*/') = f '0(i‘)|r)dl^, the density matrix takes the integral form 

P=MM=J J''<p(r)ip*(r')lr)(r'ldVdV'. (5) 

In the Appendix it is shown explicitly that, for a diagonal p, the von Neumann 
entropy then becomes 

EntvN = —Trlplnp] = — J' J ip*{r)‘ifj{r')\n[tp*{r)'ijj{r')]dVdV'. (6) 

As it must, -g^EntvN = 0. A direct calculation (see the Appendix) reveals that 

EntvN = -JJdVdV' 

\/p(r)p(r') [in ^p(r)p(r') cos 'S(r)-S(r ) ) ^ S(r)-S(r ) ^ ( S(r)-S(r ) )] ^ 

(7) 

Eq. © illustrates that the conserved von Neumann entropy is fundamentally 
different from the Boltzmann entropy of Eq. (4). Svn is a measure of the spatial 
correlation between different properties within the fluid. Furthermore, it requires 
explicit information on the velocity potential S and not just on the fluid density 
p. To the best of our knowledge it is not related to any classical measure of fluid 
entropy, thus its conservation is applied uniquely to the Madelung quantum fluid. 
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3 The Boltzmann Entropy production and compressibility effects 


Within the definition of (4), the continuity equation (2) yields 


d s 

— + V • (sUa) = pV • Ua, 


( 8 ) 


and if all fluxes vanish at the domain boundaries, we immediately obtain that 

-EntB =kB J pV- UadV = kB{V- Ua) = kB (v^5) . (9) 


dE 


Hence, the total entropy production is equal to the expectation value of the di¬ 
vergence, i.e., the Boltzmann entropy grows through expansion of the fluid and 
may decay through compression. For the Madelung fluid to be compressible, the 
quantum action S must have a “source” in the sense that it has to satisfy some 
Poisson equation in the form of V^S' ^ 0. From the wave function perspective, 
compressibility results from superposition. For a single plane wave solution of the 
form tp = Ua = VS = ^k. Thus, the advective velocity 

is simply proportional to the wavenumber k and V • Ua = 0. However, when two 
plane waves or more interfere, V • Ua 0 in general. 

In the next two simple examples we show how a superposition of plane waves 
triggers entropy growth. We will furthermore see how compressibility may lead to 
a reversible process. In the two cases we consider, the dynamics of a ID Gaussian 
density solution has been derived by m, 


p{x,t) 



( 10 ) 


Substituting (10) in the continuity equation (2) yields 


Ua{x,t) 


= X 


diner 


V • Uq 


diner 


( 11 ) 


(i) First we consider the case of a free particle. Substituting (10) and (11) in 
the ID version of (3) for U = 0 yields 


’ dt^ 


2ma 


2 2 , 
er = erp -h 


\ 2mao J 


with the Boltzmann entropy 


EntB = kB In (er\/27re) = EntBo + ^ In 


1 + 


ht 


^2maQ 


( 12 ) 


(13) 


so that 


1 dEntB „ diner 

kB dt dt 



> 0 . 


(14) 


These results illustrate how the superposition of plane waves in the Gaussian 
wave packet of Eq. (nni) influences the compressibility of the Madelung fluid and 
how it gives rise to an increase of entropy. For a free particle, a ID plane wave 
solution of (1) has the form oiE = , hence the dynamic Gaussian 
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solution of (10) is a continuous superposition of plane waves whose amplitude 
A{k) oc , as can be verified from the Fourier transform of (10) at time 

t = 0. 

(ii) As a second example we solve (10) in the presence of the harmonic potential 
U = ^{ujqx)^ to obtain from (3): 


2, 2 2n 

-cr 


2 

(To = 


2mujo 


(15) 


which can be solved numerically. The well known stationary ground state solution 

_ 

(e.g. [13]) in which a = ao, is a special case of (15) where po = i ^md 

the Boltzmann entropy is constant. 


EntBo 


-ks 


Po In podx = fcs In ((To's/27re). 


(16) 


As pointed out by [7], and is evident from (11), the only possible solution for the 
velocity in this case is Ua = 0, hence the Madelung fluid is obviously incompressible 
in the ground state. Consider, however, a small deviation from the ground state: 
a = (7Q + e(t), where |e|/(To << 1. Eq. (15) then yields ^ = —2a;Qe + O(e^), so 
that e(t) = €0 cos{\/2ojot) for 0(e). Hence, 


TP , TP , , 7 / \ 1 dEnts „ 1 5e %/2wo . , /k 

EntB = EntBo + ^b )—), - tt :— = V • Ua = —tt: =-sm(V2wot), 

(To kB at (To at (To 

(17) 

implying a reversible sinusoidal variation of the Boltzmann entropy. 

For completeness we note that the action S, associated with the wave function 
of (10), can be found explicitly (up to some constant) using Ua = x-^\na = 
-^S so that S = ^ + f{t). The time-dependent function /(t), must satisfy 

the Hamilton-Jacobi equation (or the time-dependent Bernoulli equation in the 
fluid dynamics language 0), which is the imaginary part of (1) from which (3) is 
derived. For the ID version this equation becomes 


^ 1 
dt ^ 2 



2 

Q U = 0, 


(18) 


yielding /(i) = -(^)2 . 


4 The role of the Bohm potential in classical diffusion 

We wish to compare the entropy dynamics of the Madelung fluid to the irreversible 
processes obtained in classical diffusion. We consider the standard case where the 
diffusive velocity Ud satisfies Fick’s first law: 

Ud = -DVlnp, (19) 

where D is the diffusion coefficient (assumed constant for simplicity). Changes 
in density result then from diffusive fluxes (rather than advective fluxes in the 
hydrodynamic continuity equation of (2)) as stated by Fick’s second law: 

|^ = -V-(pUd) = DV^p. 


( 20 ) 
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Equation (8) still holds when Ua is replaced by u^, but its domain integration 
yields now: 


— {EntB)D = ks 


pV • \iddV = ks {y ■ Ud) = fcs-D ^(Vlnp)^^ > 0, 


( 21 ) 


where Fi = f dV = ((Vlnp)^) is the Fisher information and {EntB)D rep¬ 
resents the Boltzmann entropy undergoing a diffusion process. Hence, this well 
known relation (e.g. m, and other information theory contexts) suggests that en¬ 
tropy increases in diffusive processes through expansion, as in the hydrodynamic 
case, but in contrast with the latter, the entropy increases irreversibly with time 
as long as density gradients exist (in agreement with the second law of thermody¬ 
namics) . 

Equation (20) can be transformed into a momentum like equation when defin¬ 
ing a “material diffusive derivative” as ^ V. Then it is straightforward 

to show that (20) can be translated to 


-^Ud = VQd, ( 22 ) 

where Qd = —2D^ may be denoted as a diffusive Bohm potential. Hence, the 
positive gradient (as opposed to the negative sign in (3)) of the diffusive Bohm 
potential acts as a “force” to accelerate the diffusion. This appearance of the 
quantum Bohm potential in a classical process is intriguing. 

The role of the gradient of the Bohm potential as a diffusive force becomes 
more transparent when returning to the dynamic ID Gaussian example. The well 
known solution to the ID version of (20) (e.g., [Tl]) is given by 


p{x,t) = 




e ^ a = 2Dt, Ud{x,t) = —D 


d\np _ Slncr 


dx 


= X- 


dt 


Hence 




X 


2t' 

(23) 


(24) 


(note that for large t both (14) and (24) experience asymptotic entropy growth 
which is proportional to t~^). Since Ud = ^ = x/{2t) a fluid element located at 

Xq at time to will be drifted at time t to a; = (random walk). Therefore, 

Ud = xl{2t) = ^acceleration of the fluid element is 


d _ j xoy/t\ ^ _ X 
) 4t2 =“4t2’ 


(25) 


where -^Ud = ^Ud = + xid^g^- Thus, in this example the diffusion rate is 

being decelerated by the gradient of the diffusive Bohm potential 


Dd _ dQd _ d f X \'2 _ X 

Dt dx dx / 41 ^ 


( 26 ) 
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5 Representation of diffusion in the Madelung fluid 

The authors of [slleirfl defined the complex velocity derived from the momentum 
operator —as 

V = [—i—V Inif/] = Vr-+ fvi, (27) 

m 

so that -Vr = Ua = VS is the advective velocity and Vi = — ^V(lnp). One may 
interpret = u^, as in (19), suggesting the relation between the Planck constant 
and the diffusion coefficient to be ^ = mD, so that Qd becomes identical to the 
quantum Bohm potential. For the simplest case of the Einstein relations in an 
ideal gas, mD = (1 • p)/3, where |1| is the molecular mean free path and |p| is 
the magnitude averaged thermal (random walk) molecular momentum in-between 
collisions. Isotropy results in the relation ^ = IxPx, and indeed lx and px are the 
basic scales obtained from statistical mechanics for deriving the kinetic theory of 
gases. In other words, one cannot resolve the ideal gas dynamics within length 
scales smaller than lx or for momenta smaller than px- As pointed out by |15l 
[T6| , it is intriguing that these two fundamental scales form the canonical variables 
which set the exact limiting case of the Heisenberg uncertainty principle. 
Incorporating (2) with (19) we obtain the Fokker-Planck equation 

f + V. [p(u.-„,)] = (A) W (28) 

which corresponds to the entropy equation 

-h V • [(s - p)Ua] = pUa • Ud, (29) 

yielding 

= ’^B{ua-nd). (30) 

Thus stating that positive entropy production in the Madelung fluid occurs when 
the advective and the diffusive velocities are positively correlated within the fluid 
domain. 


6 Summary and Discussion 

Having an intuitive classical interpretation, the Madelung formulation is natural 
for discussing quantum dynamics and thermodynamics. In Ref. [7], two of us high¬ 
lighted the role of compressibility in linking hydrodynamical and thermodynamical 
processes in the Madelung fluid. 

In the current work, when examining entropy production in this framework, 
we And that it is proportional to the expectation value of V • Ua. Thus, the ex¬ 
pansion of the fluid is equivalent to an entropy increase. Furthermore, we note 
that interference between plane waves leads to an increase of entropy, but it also 
renders the Madelung fluid compressible, allowing it to expand. Putting all of 
the pieces together, we find that a simple link exists between the interference of 
wave functions, compressibility of the Madelung fluid, and (Boltzmann) entropy 
production. 
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We illustrated how these concepts come to life in several specific examples. We 
underscored the difference between the Boltzmann and von Neumann entropies of 
the Madelung fluid (and derived simple new form for the latter). 

The Madelung fluid expansion suggests an analogy with diffusive processes 
such as a free Joule expansion. Ref. [7] related the imaginary part of the quantum 
velocity to the thermal fluctuations of the Madelung fluid. Following the current 
analysis, this imaginary part can indeed be interpreted, as a diffusive drift velocity, 
since its flux is propositional to the minus sign of the density gradient (as in 
Pick’s first law). Furthermore, this analogy suggests that the diffusion coefficient 
is /i/ 2 m, and for the simplest model of diffusion in ideal gas, the relation between 
the diffusion coefficient, the mean free path, and Brownian momentum provides 
the lower bound of Heisenberg’s uncertainty relation. 

The entropy production can be expressed as well in terms of the correlation 
expectation value between the advective and diffusive velocities. Moreover, for 
classical diffusive processes, it was shown that the gradient of the Bohm potential 
acts as a force to accelerate the diffusion. This appearance of the quantum poten¬ 
tial in a classical mechanism may shed light on its role in the Schrodinger equation. 
It is important to remember however, that in the quantum case the diffusion equa¬ 
tion, governed by the Pick’s second law, cannot be extracted from the Schrodinger 
equation, i.e., ^ 7 ^ {h/2m)V^p, but rather ^ = — V - (pUa). Therefore, processes 
associated with density variation, such as entropy growth, are due to expansion by 
the hydrodynamical advective velocity and not by diffusive irreversible processes. 
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Appendix 

Let our system be described by a continuous wavefunction lip) = f dx'ip{x)\x) 
(to simplify the notation we shall use in the appendix a 1-dimensional system, 
although the results are completely general). Then the pure density matrix would 
be: 

p = |if'){!P'| = J j dxdx'i'*{x)i'{x')\x){x'\, (31) 

with entries 'P* {x)^{x'). Therefore, in a basis where p is diagonal 
plnp = 

= 1111 dxdx'dydy'P*{x)<P{x')ln[P*{y)<I'{y')]\x){x'\y){y'\ = 

— fill dxdx'dydy'P*{x)P{x')ln[P*{y)^{y')]5{x' — y)\x){y'\ = ' ' 

= 111 dxdydy''I'*{x)'P{y)ln[<P*{y)<P{y')]\x){y'\ = 

To evaluate the von Neumann entropy © we shall now multiply the previous 
equation by {z\ from left, by \z) from right and then integrate over z do find: 

EntvN = - j j J J dxdydy'dz'P*{x)'P{y)ln[P*{y)'I'{y')\{z\x){y'\z) = 

= — f f f f dxdydy'dzP*(x)P(y)ln[>P*(y)>P(y')]S(x — z)5{y' — z) = 

= - f f dydzP*(z)^(y)ln[P*(y)P(z)] = 

= —ff dxdx'<P*(x)<P(x')ln[<P*(x')<plx)], 


(33) 
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where we only changed variables in the last line. 

The von Nenmann entropy of an isolated system is constant nnder unitary 
time evolution (this follows, e.g., from the cyclicity of the trace). We will show 
that explicitly: 


_aEr^ = r { dxdx'{dt^*(xW(x')lnW*{x')'P{x)U 
{x)dti'{x')ln[I'* {x')’P{x)] + 




(34) 


We now use the Schrddinger equation for evaluating all the time derivatives. Next, 
we differentiate twice by parts the first term. It cancels the last term while leaving 
an additional —i'I'*{x)'P{x') (up to a multiplicative constant). Similarly, 

when differentiating twice by parts the second term, it cancels the third term, 
leaving an additional Therefore, 


dEntvN 

dt 


i f f dxdx'P*{x)'P{x') 




Int. 


(35) 


The von Neumann entropy and its derivatives are real, hence Int = Int*. However, 
if we perform a change of variables in Int* such that x ^ x' and x' x, we find 
out that Int* = —Int. Therefore, Int = 0 as required. 
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